Observation of an enhancement in e+e — t- T(15')7r+7r , T(2S')7r+7r , and T(3S')7r+7r 

production around y/s = 10.89 GeV at Belle 
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We measure the production cross sections for e^e~ — >■ T(15)7r+7r~, T(25')7r"'"7r~ , and T(3S)7r"'"7r~ 
as a function of ^/s between 10.83 GeV and 11.02 GeV. The data consists of 8.1 fb~^ collected with 
the Belle detector at the KEKB e'^e~ collider. We observe enhanced production in all three final 
states that does not conform well with the conventional T(10860) lineshape. A fit using a Breit- 
Wigner resonance shape yields a peak mass of [I0888.4t^;^(stat) ± 1.2(syst)] MeV/c^ and a width 
of [30.7t?;g(stat) ±3.1(syst)] MeV/c^. 

PACS numbers: 13.25.Gv, 14.40.Pq 



A host of new charmonium-like mesons have been 
discovered recently that do not seem to fit into the 
conventional cc spectrum. Possible interpretations for 
these exotic states include multiquark states, mesonic- 
molecules (cq-cq, where q represents a u-, d- or s-quark), 
or ccg "hybrids" (where is a gluon). The narrow 
X(3872) m, and various vector mesons, in particular, the 
broad y(4260) (i-Q resonance, were revealed by their di- 
pion transitions to J/V' (and i}}'). By analogy with the 
y (4260) discovery, it was suggested that a hidden-beauty 
counterpart, denoted Yf, |5;], could be searched for in ra- 
diative return events from e^e~ at center-of-mass (CM) 
energy on the T(10860) peak, or by an energy scan above 
it. 

Analyzing a data sample recorded by the Belle detec- 
tor at a single CM energy ^/s near the T(10860) res- 
onance peak, no unusual structure was observed below 
the peak. Instead, anomalously large T(lS')7r+7r^ and 
T(2iS')7r+7r~ production rates were observed p. If these 
signals are attributed entirely to dipion transitions from 
the T(10860) resonance, the corresponding partial widths 
would be two orders of magnitude larger than those for 
corresponding transitions from the T(2S'), T(3S'), and 
T(4S') states. Possible explanations include a new non- 
perturbative approach [3] to the decay widths of dip- 
ion transitions of heavy quarkonia, the presence of fi- 
nal state interactions [1| , or the existence of a tetraquark 
state [l3| • A measurement of the energy dependence of 
the cross sections for e+e" — > T(n5)7r+7r~ {n — 1,2,3) 
in and above the T(10860) energy region will provide 
more information for exploring these models. 

Here we report the observation of the production of 
e+e- T(lS')7r+7r-, T(2S')7r+7r-, and T(35)7r+7r- at 
^/s ~ 10.83, 10.87, 10.88, 10.90, 10.93, 10.96, and 11.02 
GeV, while the previous data collected at ^/s ~ 10.87 



GeV was already reported in Ref. 

The new 8.1 fb~^ data sample was also collected with 
the Belle detector at the KEKB e+e~ colhder [Ilj. In 
addition, nine scan points with an integrated luminosity 
of ~ 30 pb~^ each, collected in the range of -y/s = 10.80- 
11.02 GeV, are used to obtain the hadronic line shape of 
the T(10860). 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-fiight scintillation counters 
(TOP), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect if£ mesons and to identify muons (KLM). The 
detector is described in detail elsewhere [l2l |. 

Events with four well-reconstructed charged tracks and 
zero net charge are selected to study T(niS')7r"'"7r~ pro- 
duction. A final state that is consistent with a T{nS) 
candidate and two charged pions is reconstructed. A 
T{nS) candidate is formed from two muons with op- 
posite charge, where the muon candidates are required 
to have associated hits in the KLM detector that agree 
with the extrapolated trajectory of a charged track pro- 
vided by the drift chamber. The other two charged 
tracks must have a low likelihood of being electrons 
(the likelihood is calculated based on the ratio of ECL 
shower energy to the track momentum, dE/dx from the 
CDC, and the ACC response); they are then treated as 
pion candidates. This suppresses the background from 
e+e^ — >■ ^"'"/i~7 — > /i+/i~e'^e" with a photon conver- 
sion. The cosine of the opening angle between the tt"*" 
and TT^ momenta in the laboratory frame is required to 
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TABLE I: CM energy (\/s), integrated luminosity (£), signal yield (Ns), reconstruction efficiency, and measured cross section 
(a) for e+e^ — )■ T(15')7r+7r^, T(25')7r^7r~ , and T(35')7r+7r~. Due to a negative yield, an upper limit at 90% confidence level 
for a[e+e- T(3S)7r+7r-] at = 10.9555 GeV is given. 



v^(GcV) £(fb-i) 



3.3 
+ 7.2 
6.5 
+5.8 
5.1 



e+e- -s. T(15)7r+7r- 
Efr.(%) a(pb) 



0.56lo^g ± 0.06 
2.14l0 j6 ± 0.15 
l.6S+°il ± 0.13 
0.89t°;27 ± 0.08 
0.37t° i| ± 0.04 
0-53lo.23 ± 0-05 



e+e- -s> T(2S)7r+7r- 
Afs Eff.(%) a{ph) 



+b.b 
4.9 
+9.0 
-8.3 
+7.4 
6.7 
+3.8 
3.1 
+2.0 
1.3 
+3.1 
2.4 



2.05lj;'*5 ± 0-24 
llf6l±0.59 
4.53t!5-^7 ± 0.51 
1.19l»j^±0.16 
0.29lo ig ± 0.05 

^+o:5i 



1.1 

+4.3 
3.7 
+3.7 
3.1 
+2.2 
1.5 
+2.5 
-3.0 
+2.6 
1.9 



Efr.(%) 



T(3S)7r+7r- 
a{ph) 



0.23lo;i4 ± 0.03 
1.47155'*? =•= 0-18 
1.26tEj;« ± 0-15 
0.4ll55;^j ± 0.05 
-0.28t°-47 ±0.03 < 0.20 

-,+0.42 . 



1.73 
1.89 
1.46 
1.18 
0.99 
0.88 



10.6: 
43. 4 j 
26.2j 



3.9 
4.9 



3.3 
+2.6 
1.9 
+2.8 
2.1 



43.8 
43.1 
43.2 
42.6 
42.5 
42.0 



24.0 

68 

45.4 
9.7; 
2.0 
5.5 



34.9 
35.4 
35.6 
35.9 
36.4 
36.0 



5.31^ 



1.8 
14.9 
10.3 

2.9 
-1 

4.3 



20.5 
24.5 
25.7 
27.5 
29.4 
32.7 



0.69" 



19 



+4.U 
3.3 



0.30 



± 0.08 



+U.bb 



21.74 



325 



37.4 1.61 ± 0.10 ± 0.12 



186 ± 15 18.9 2.35 ±0.19 ±0.32 



10.5 



1.5 
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± 0.19 
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FIG. 1: The distributions of AM- [^-Mt(„s)] = 1,2,3) for (a-f) T(15)7r+7r-, (g-1) T(2S)7r+7r-, and (m-r) T(35')7r+7r- 
events with the fit results superimposed. The six columns of plots represent the data samples collected at different CM energies. 
The dashed curves show the background components in the fits. 



be less than 0.95. The four-track invariant mass must 
satisfy \M{fi+fi-T:+T:-) - V^l < 150 MeV/c^. The trig- 
ger efficiency for four-track events satisfying these criteria 
is very close to 100%. 

The kinematic variable AM, defined by the differ- 
ence between M{^^ ^~Tr'^Tr~) and M{^'^ii~), is used 
to identify the signal candidates. Sharp signal peaks 
are expected at AAf = \fs — M-x{nS)- The candidate 
events are separated into three distinct regions defined 
by |AAf - [yi - Mx(„s)]| < 150 MeV/c^ for n = 1, 2, 
and 3, and signal yields are extracted from an unbinned 
extended maximum likelihood (ML) fit to the AAf dis- 
tribution within each region. The likelihood function for 
each fit is defined as 

L{Ns,Nt) = —^^ Y[[Ns-Ps{AM,)+Nb-Pb{AMi)] , 

i=l 

where Ng {Nb) denotes the yield for signal (background), 
and Ps {Pb) is the signal (background) probability den- 
sity function (PDF). The signal is modelled by a sum of 
two Gaussians while the background is approximated by 



a linear function. The Gaussians parameterized for the 
signal PDF are fixed from the Monte Carlo (MC) simu- 
lation at each energy point. We fit 18 AM distributions 
(shown in Fig. [1] with the fit results superimposed) si- 
multaneously with common corrections to the mean and 
width of the signal Gaussians. 

The measured signal yields, reconstruction efficiencies, 
integrated luminosity, and the production cross sections, 
as well as the results from the previous publication Q 
for the data sample collected at i/s = 10.867 GeV, are 
summarized in Table ID The efRciencies for T(35)7r"'"7r~ 
are much improved compared to Ref. @ as the inefficient 
selection criterion ^max < 175° has been removed, where 
^'max is the maximum opening angle between any pair of 
charged tracks in the CM frame. 

For the cross section measurements, systematic uncer- 
tainties are dominated by the T{nS) — > fi^ branch- 
ing fractions, reconstruction efficiencies, and PDF pa- 
rameterization for the fits. Uncertainties of 2.0%, 8.8%, 
and 9.6% for the T(15'), T(2S'), and T{3S) fi- 
branching fractions are included, respectively. For the 
T(lS')7r+7r~ and T(2S')7r+7r~ modes, the reconstruction 
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efficiencies are obtained from MC simulations using tlie 
observed M(7r"'"7r~) and cos^Hei (the angle between the 
7T~ and T(10860) momenta in the tt+tt^ rest frame) 
distributions in our previous publication as inputs 
The uncertainties associated with these distributions give 
rise to 2.7%-4.5% and 1.9%-4.2% uncertainties for the 
T(lS')7r+7r~ and T(2S')7r+7r^ efficiencies, respectively. 
The ranges on the uncertainty arise from the CM en- 
ergy dependence of the tt+tt^ system. We use the model 
of Ref. [l^l as well as a phase space model as inputs 
for T(35')7r"'"7r~ measurements; the differences in accep- 
tance are included as systematic uncertainties. The un- 
certainties from the PDF parameterization are estimated 
either by replacing the signal PDF with a sum of three 
Gaussians, or by replacing the background PDF with a 
second-order polynomial. The differences between these 
alternative fits and the nominal results are taken as the 
systematic uncertainties. Other uncertainties include: 
tracking efficiency (1% per charged track), muon iden- 
tification (0.5% per muon candidate), electron rejection 
for the charged pions (0.1-0.2% per pion), trigger efficien- 
cies (0.1-5.2%), and integrated luminosity (1.4%). The 
uncertainties from all sources are added in quadrature. 
The total systematic uncertainties are 7%-ll%, 11%- 
16%, and 12%-14% for the T(15')7r+7r-, T(25)7r+7r-, 
and T(3S')7r^7r~ channels, respectively. 

In order to extract the resonance shape, we perform a 
fit to the measured production cross sections, includ- 
ing the one estimated at -^s = 10.867 GeV, using the 



model crr{nS)i77r/<7° 



^T(nS)7r7r 



\Ro + e"l'BW{^l,T)\ 



where = Aira^ /is is the leading-order e+e 
cross section and BW{fi, F) is the Breit-Wigner function 
l/[(s — ^^) + z/iF]. The normalizations for T(lS')7r"'"7r~, 
T(2S')7r+7r^, and T(3S')7r+7r~ {Ax(nS)TrTr), as well as the 
amplitude of the flat component Rq, the mean fj,, width F, 
and the complex phase of the parent resonance are free 
parameters in the fit. Because of the low statistics, com- 
mon resonance parameters are introduced for the three 
different final states. Results of the fits, shown as the 
smooth curves in Fig. [51 are summarized in Table [ill 
The fit quality is — 24.6 for 14 degrees of freedom. 
An alternative fit without the last data point collected at 
Vs ^ 11.02 GcV yields a similar result, ^ = 10889.0t^:^ 
MeV/c^ F = Sltll MeV/c2, and 21.3 for 11 de- 

grees of freedom. Systematic uncertainties associated 
with the cross section measurements are propagated to 
the resonance shapes. The fits are repeated, and the 
variations on the shape parameters are included as the 
systematic uncertainties. In addition to the uncertainties 
on the cross sections, the beam energy around T(10860) 
is measured by Mx(„5) -I- AM in the T(nS')7r"'"7r~ events, 
and an uncertainty of ±1 MeV is included. For the scan 
data, a common energy shift is also obtained from the 
fit to T(nS')7r+7r~ events. The relative beam energies 
are further checked using the M{^+n~) distributions of 
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FIG. 2: The CM energy-dependent cross sections for e e~ — >■ 
T(n5')7r^7r^ (n — 1, 2, 3) processes normalized to the leading- 
order e+ e~ — >■ fi+ fj,~ cross sections. The results of the fits are 
shown as smooth curves. The vertical dashed line indicates 
the energy at which the hadronic cross section is maximal. 



/i-pair samples. 

The resonance parameters for the T( 10860) are de- 
termined using energy scan data collected at CM en- 
ergies between 10.80 and 11.02 GeV. We measure the 

<^<°-2(s)/£e6(5) is 



fraction Ri, = Ub/a^^, where ab 
the e+e^ — > bb hadronic cross section. The number of 
e+e" -> bb events with i?2 < 0.2 (7V^^<°-^) is estimated 
by subtraction of non-66 events scaled from a data set 
collected at y/s ~ 10.52 GeV, where i?2 denotes the ra- 
tio of the second to zeroth Fox- Wolfram moments [l^. 
Selection criteria for hadronic events are described in 
Ref. 15|. The acceptance for e+e~ -> 66 {eb{s)) is found 
to vary slightly from 68.1% to 70.5% over the range of 
scan energies. The line shape used to model our data is 
given by I |2 -f I Ao+Awseoe"*""'"" BW {jjioseo , Tioseo ) + 
^110206''^"°"' ^W^(A*ii020,riio2o)P; this parameteriza- 
tion is the same as that used in Ref. [IQ] . We perform a 
X^ fit to our Rb measurements as shown in Fig.[3l[a). The 
shapes for T(11020) {(j)rw2o, M11020 , and F11020) are fixed 
to the values in Ref. [16[, since our data points are not 
able to constrain the T(11020) parameters. The resulting 
shape parameters for the T(10860) are ^loseo = 2.331q 24 



TABLE II: Cross sections (a) at peak, mean (/i), width 
(F), phase {(j)), and the amplitude for the constant compo- 
nent (Ro) from the fit to the CM energy-dependent e+e^ — >■ 
T(15')7r+7r~ , T(2S')7r+7r~ , and T(3S')7r"'"7r~ cross sections. 
There are two solutions for (f) and Ro with identical ■ The 
first uncertainty is statistical, and the second is systematic. 



T(15')7r7r a at peak 
T(2S')7r7r a at peak 
T(35')7r7r a at peak 



(2.78to;34 ± 0.23) pb 
(4.82l[J:^^ ± 0.66) pb 
(1.7 11°:^^ ± 0.24) pb 



^ (I0888.4t^;^ ± 1.2) MeV/c^ 

r (30.7l?:g ±3.1) MeV/c^ 

(j> (1.97 ± 0.26 ± 0.06) or (-1.74 ± 0.11 ± 0.02) rad 

Ro (l-98to.60 ± 0-20) or {0.87toil ± 0.09) (GeV)-^ 
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FIG. 3: (a) The Rt and (b) the ratio between a[e^e^ — )■ T(nS)7r+7r^] and (T[e+e^ — >■ bb] as a function of CM energy; (c) the 
energy-dependent cross section ratios for e'^e~ — > T(nS')7r^7r~ events, the result of fits with resonant parameters from Rb or 
PDG averages are superimposed. The horizontal dotted line in (a) is the non-interfering |j4„rp contribution in the fit. The 
vertical dashed line indicates the energy at which the hadronic cross section is maximal. 



rad, /Z10860 = 10879 ± 3 MeV/c^ , and Tioseo = 46t? 
MeV/c^. These values are consistent with those obtained 
in Ref. [ll]. The quality of the fit is = 4.4 for 9 degrees 
of freedom. 

Figure ^h) shows the ratio between a[e^e^ — > 
T(nS')7r+7r^] and (j[e^e^ — > bb] as a function of CM 
energy. A fit to the T{nS)TnT cross sections with 
MT(nS)7r7r = M10860 and rx(„s)7r7r = Tioseo givcs = 
37.8 for 16 degrees of freedom as shown in Fig.[3jc). This 
corresponds to a deviation of 3.2(t from the nominal fit 
with floated /iT(nS)7rir and TY(nS}Ti-w If we perform a scan 
within the region with {^J.T(nS)^T7v - Mio86o)^/ct(/^io86o)^ + 
(rT(«s)7r7r " rio86o)^/o-(rio86o)^ < 1, a minimum devia- 
tion of 2.3(7 is found, where o'(/iio86o) and (T(rio86o) are 
statistical. If the resonant mean and width are fixed to 
the results in Ref. [l6[ or PDG values jT3|, a deviation of 
3.9(7 or 5.6a is obtained; scans within the la bound yield 
a deviation of 3.4(7 or 5.1a. 

In summary, we report the observation of enhanced 
e+e- T(lS')7r+7r-, T(2S')7r+7r-, and T(35)7r+7r- pro- 
duction at CM energies between y/s ~ 10.83 and 11.02 
GeV. The energy- dependent cross sections for e+e^ 
T(nS')7r+7r~ events are measured for the first time, and 
are found to differ from the shape of the e'^e~ — ^ bb 
cross section. A Breit-Wigner resonance shape fit yields 
a peak mass of 10888.4+2;^ (stat) ± 1.2 (syst) MeV/c^ 
and a width of 30.7l:? i^ (stat) ± 3.1 (syst) MeV/c^. A 
fit excluding the y/s ^ 11.02 GcV data point is consis- 
tent with the nominal fit, indicating no strong contribu- 
tion of T(n5)7r7r events from T(11020). The T(10860) 
shape parameters obtained from our Rb hadronic cross 
section measurements are consistent with the measure- 
ments from BaBar. The differences between the shape 
parameters from T(n5)7r7r events and from our Rf, mea- 
surements are iir(nS)-Kn — /^loseo = 9 ± 4 MeV/c^ and 
^Tins)-.^ - rio86o = "ISt^a MeV/c^. A fit to the 
e+e^ — > T(n5)7r+7r^ cross sections with our measured 
T(10860) mean and width yields a deviation of 3.2ct, 
while a deviation of 2.3(7 is obtained if we allow a max- 
imal ±1(7 drift of the T(10860) shape parameters. The 
T(nS')7r"'"7r~ partial widths were found to be much larger 
than the expectations for conventional T{5S) states. As 



an extension, energy- dependent T(nS')7r^7r~ production 
cross sections are measured; the observed structure devi- 
ates from the T(10860) shape obtained from the hadronic 
cross sections. 
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